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1.  Introduction 

Oscillations  of  an  isolated  thin  magnetic  tube  at  various  depths  of  the  convective  zone  and 
in  the  atmosphere  of  the  Sun  are  investigated  in  the  present  work  in  the  approximation  of  ideal 
magnetic  gas-dynamics  ([1]).  These  oscillations  are  the  source  of  generation  of  acoustic  waves 
whose  transition  in  weak  shock  waves  provides  anomalous  heating  of  the  solar  atmosphere. 

The  modes  of  small-amplitude  linear  oscillations  at  various  depths  of  the  convective  zone 
are  investigated.  The  loss  of  stability  of  equilibrium  of  the  magnetic  tube,  which  is  initially  at 
rest  in  the  convective  zone,  as  well  as  the  nonlinear  large-scale  oscillations  of  the  tube  and  tube 
emergence  into  the  solar  atmosphere  are  independently  investigated. 

The  results  obtained  are  fundamental  for  research  of  sunspot  generation  and  physics  of  so¬ 
lar  activity;  they  are  important  for  analysis  of  generation  and  time  evolution  of  global  magnetic 
structures,  solar  wind,  and  some  other  astrophysical  phenomena  [2]. 

2.  Mathematical  statement  of  the  problem 

In  this  work,  radial  oscillations  of  a  thin  magnetic  tube  located  in  the  equatorial  plane  of 
the  Sun  (Fig.  la)  are  investigated.  The  zeroth  harmonic  of  oscillations  is  studied,  at  which  only 
the  magnetic  tube  radius  varies  (Fig.  lb). 

The  full  set  of  equations  of  ideal  magnetic  gas  dynamics  [1,  3]  is  shown  below: 

2iirpE  =  M0  =  const ,  (1) 


Figure  la.  Location  of  a  thin  magnetic  tube  in 
the  equatorial  plane  of  the  Sun. 


Figure  lb.  System  of  coordinates  for  calculations. 
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Equation  (1)  formulates  the  law  of  mass  conservation  in  the  magnetic  tube.  Equations  (2) 
and  (3)  describe  the  tube  motion  under  the  influence  of  the  archimedean  force  and  tension  of 
magnetic  field  lines.  The  velocity  of  the  tube  is  supposed  to  be  essentially  subsonic,  and  the 
influence  of  the  ambient  medium  on  the  tube  motion  is  expressed  in  terms  of  additional  mass  in 
(3)  [4].  Equation  (4)  is  the  adiabatic  energy  law.  Equation  (5)  is  the  law  of  conservation  of 
magnetic  flux  over  the  tube  cross  section.  All  parameters  depend  only  on  the  radius  of  the  tube 
curvature  -  r.  Equation  (6)  determines  the  equality  of  the  ambient  gas  pressure  pe  to  the  internal 
pressure  in  the  magnetic  tube,  which  consists  of  the  gas  pressure  and  the  pressure  of  the  mag¬ 
netic  field.  The  dependences  of  pe(r),  g(r),  pe(r)  on  r  are  determined  by  the  chosen  model  of  the 
solar  convective  zone  [5],  and  they  are  necessary  to  close  the  system  of  equations. 

System  (l)-(7)  is  dimensionless,  and  the  units  of  plasma  parameters  in  the  tube  are  shown 
below: 

Do]  =  1010  cm  [go]  =  1.327- 106  cm/s2  [p0]  =  105  dyn/cm3  [r0]=104K 

[po\  =  10“6  g/cm3  [t0]  =  (r0/g0)1/2  [u0]  =  g0-t0 

3.  Liner  oscillations  of  the  magnetic  tube  at  various  depths  of 
the  convective  zone  of  the  Sun.  Relaxation  zone  and  zone  of  solar  Dynamo 

The  equilibrium  position  of  the  magnetic  tube  can  be  achieved  by  imposing  the  following 
limitation  on  gas  density  and  magnetic  field  strength: 

1  H2 

(P-P  e)g  =  ~ - •  (8) 

4ji  r 

In  equilibrium,  the  force  of  magnetic  field  tension  counterbalances  the  archimedean  buoy¬ 
ancy  force.  In  equilibrium,  the  magnetic  field  strength  corresponding  to  (6),  (7),  (8)  and  (3) 
uniquely  determines  the  plasma  density,  pressure  and  temperature. 

The  maximum  possible  value  of  magnetic  field  intensity  in  the  tube  is  determined  from 
Eq.  (6): 

#max  =  V871^  •  (9) 

If  the  depth  of  tube  location  is  less  than  40  km  below  the  photosphere,  then  the  tube  in 
equilibrium  is  unstable  for  all  possible  values  of  the  field  strength  (9).  The  tube  located  in  the 
range  of  depths  from  40  km  to  140- 103  km  below  the  photosphere  is  in  stable  equilibrium  if  the 
magnetic  field  intensity  is  greater  than  the  critical  value  Hcr(r)  (Fig.  2).  The  tube  equilibrium  in 
that  range  of  depths  is  unstable  if  the  magnetic  field  intensity  is  lower  than  the  critical  value. 


(2) 

(3) 
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H,  10  6Gs 


Figure  2.  Critical  intensity  of  the  magnetic  field 
versus  the  depth  in  the  relaxation  zone. 


Hcr,  10  7Gs 


-210  -205  -200  -195  -190 

Figure  3.  Critical  intensity  of  the  magnetic  field 
versus  the  depth  in  high  layers  of  the  Dynamo 
zone 


The  second  area  where  the  tube  can  be  in  stable  equilibrium  is  situated  lower  than  190-103 
km  below  the  photosphere  (Fig.  3).  In  contrast  to  the  previous  domain,  the  tube  here  is  in  steady 
equilibrium  if  the  magnetic  field  intensity  is  lower  than  the  critical  value  Hc,.(r)  depending  on 
depth  (see  Fig.  3).  If  the  magnetic  field  intensity  is  greater  than  the  critical  value,  equilibrium  of 
the  magnetic  tube  becomes  unstable,  and  the  tube,  deviating  from  its  initial  position  due  to  a 
small  perturbation,  emerges  upward  with  an  increasing  velocity  due  to  convective-like  instabil¬ 
ity  development. 

Lower  than  209T03  km  below  the  photosphere,  the  critical  values  of  the  magnetic  field  in¬ 
tensity  exceed  the  maximum  possible  value  determined  by  (9).  Below  this  depth,  there  are  ideal 
conditions  for  storage  and  generation  of  magnetic  energy.  The  region  below  209T03  km  is 
termed  as  "the  zone  of  action  of  the  solar  Dynamo"  or,  simpler,  "the  Dynamo  zone  "  [6]. 


4.  Emergence  of  magnetic  fields  from  the  Dynamo  zone 

The  range  of  depths  (190T03  -s-  209T03  km)  where  the  generated  magnetic  fields  lose  sta¬ 
bility  and  float  up  to  the  photosphere  if  their  intensity  exceeds  the  critical  level  practically 
forms  the  majority  of  all  nonstationary  active  solar  processes  [7]. 

The  processes  of  the  loss  of  stability  and  emergence  of  magnetic  fields  in  the  upper  layers 
of  the  convective  zone  are  nonlinear;  therefore,  their  investigation  requires  numerical  methods 
to  solve  system  (l)-(7).  At  the  initial  time,  the  condition  of  equilibrium  is  set,  the  tube  is  as¬ 
sumed  to  be  at  rest;  thus,  only  the  relation  between  the  field  intensity  and  the  critical  value  Hcr 
at  the  same  depth  determines  the  subsequent  evolution  of  the  tube. 

If  the  initial  magnetic  field  intensity  of  the  tube  is  greater  than  the  critical  one,  the  tube 
equilibrium  loses  stability  and  the  tube  is  floating  up  with  an  increasing  velocity  in  the  convec¬ 
tive  instability  regime.  Results  of  numerical  modeling  of  tube  emergence  out  of  the  Dynamo 
zone  from  various  initial  depths  are  shown  in  Fig.  4. 

The  temperature  and  density  of  the  tube  plasma  decrease  due  to  adiabatic  extension  of  the 
emerging  tube.  There  exists  a  special  depth  where  the  tube  acceleration  changes  its  sign  and  the 
tube  begins  to  break,  which  depends  on  the  initial  parameters  of  the  tube.  The  total  buoyancy 
force  is  equal  to  zero  at  the  point  where  the  force  changes  its  sign.  This  point  plays  the  role  of 
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Figure  4.  Tube  emergence  out  of  the  Dynamo  Figure  5.  New  position  of  the  equilibrium  point 

zone  from  various  initial  depths:  1  -  1 90- 1 03  km;  as  a  function  of  initial  positions  within  the  Dy- 

2  -  200-103km;  3  -  20  MO3  km  namo  zone 


new  equilibrium  during  emerging  of  the  tube.  The  velocity  of  emergence  is  maximal  at  this 
point.  As  a  result  of  tube  breaking,  the  velocity  of  emergence  changes  its  sign,  and  the  tube 
passes  through  the  point  of  maximum  altitude  and  begins  to  move  in  the  opposite  direction. 
Since  dissipative  processes  are  omitted,  the  tube  returns  to  the  initial  position  (Fig.  4). 

The  initial  position  within  the  Dynamo  zone  (Fig.  5)  unequivocally  determines  the  new 
position  of  the  equilibrium  point  due  to  imposing  of  the  stability  loss  condition:  H(t=0)  = 
Hcl.(r).  The  new  equilibrium  positions  are  grouped  in  a  narrow  layer  (23T03km  h-  38-103km). 
The  velocities  of  emergence  of  the  magnetic  tube  at  these  depths  are  maximal,  and  the  velocity 
can  reach  up  to  40  km/s  (Fig.  4).  The  velocity  of  emergence  of  the  magnetic  tube  remains  sub¬ 
sonic  for  all  variants  considered.  Near  the  photosphere,  the  velocities  of  the  magnetic  tube, 


v,km/s  v/Cg 
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Figure  6.  Velocity  of  emergence  of  the  magnetic  Figure  7.  Velocities  of  the  magnetic  tube,  meas- 

tube  versus  time.  ured  in  Mach  numbers,  versus  time. 

The  initial  depths  of  emergence:  1  -  190-103km;  2  -  200-103km;  3  -  201 -103  km 
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-210  -205  -200  -195  -190  -204  -200  -196  -192  -188 

Figure  8.  Full  period  of  oscillations  versus  the  Figure  9.  Direct  transition  of  magnetic  fields 

initial  depth  of  emergence.  from  the  Dynamo  zone  into  the  solar  atmosphere 

versus  the  initial  depth  of  emergence. 

measured  in  Mach  numbers,  increase  nonlinearly.  The  magnetic  tube  generates  acoustic  waves, 
which  penetrate  through  the  photosphere  in  the  solar  atmosphere  and  transform  into  weak  shock 
waves  in  the  solar  chromosphere  [8], 

The  motion  between  the  start  position  and  the  equilibrium  point  is  characterised  by  a  low 
velocity  relative  to  rapid  acceleration  and  by  changing  of  the  motion  direction  near  the  point  of 
equilibrium  (Figs.  6  and  7).  The  oscillations  have  the  character  of  separated  overshoots.  The 
duration  of  overshoot  motion  is  obviously  less  than  the  period  of  oscillations  of  the  magnetic 
tube,  and  its  characteristic  scale  is  60  minutes  (1  hour)  instead  of  the  full  period  of  oscillations 
whose  scale  is  of  the  order  of  several  days  (Fig.  8).  The  depth  of  the  maximal  ascent,  measured 
from  the  photosphere,  increases  with  decreasing  initial  depth  within  the  Dynamo  zone.  Below 
the  depth  of  203T03km,  a  direct  transition  of  magnetic  fields  from  the  Dynamo  zone  into  the 


Vph ,  km/s  v,  km/s 
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Figure  10.  Velocity  of  emergence  of  magnetic  Figure  1 1.  Magnetic  tube  velocity  distributions 
fields  at  the  photosphere  level  versus  the  initial  over  the  altitude  within  the  limits  of  the  solar 
depth  of  emergence  within  the  Dynamo  zone.  chromosphere. 
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solar  atmosphere  occurs  (Fig.  9)  [7].  The  velocity  of  emergence  of  magnetic  fields  at  the  photo¬ 
sphere  level  increases  nonlinearly  lip  to  hypersonic  values  (M  ~  4,  Fig.  10)  with  increasing  ini¬ 
tial  depth  within  the  Dynamo  zone.  Figure  11  shows  the  magnetic  tube  velocity  distributions 
over  the  altitude  within  the  limits  of  the  solar  chromosphere.  If  the  Mach  number  of  the  tube 
exceeds  2,  then  the  velocity  of  emergence  within  the  solar  photosphere  depends  weakly  on  the 
altitude.  This  circumstance  allows  us  to  explain  the  formation  of  coronal  jets  (transients)  at 
large  altitudes  within  the  solar  atmosphere  [9,  10]. 

REFERENCES 

1.  Alekseenko  S.V.,  Mezentsev  A.V.,  Romanov  V.A.,  Romanov  D.V.,  Romanov  K.V.  Stabilization  of  emerging 
magnetic  fields  in  the  upper  layers  of  the  Solar  convection  zone  //  Russ.  J.  Eng.  Thermophys.  1998.  Vol.  8.  N°  1-4. 
P.109-119. 

2.  Priest  E.R.  Solar  Magnetic  Hydrodynamic.  Moscow:  Mir,  1975.  589  p. 

3.  Spruit  H.C.,  Zweibel  E.G.  Convective  instability  of  thin  flux  tubes  //  Solar  Phys.  1979.  Vol.  62.  P.  15-22. 

4.  Landau  L.D.,  Lifshitch  E.M.  Hydrodynamic.  Moscow:  Nauka.  1986.  736  p. 

5.  Christensen-Dalsgaard  J.,  Dappen  W.,  etc.  The  current  state  of  Solar  modeling  //  Science.  1996.  Vol.  272. 

6.  Romanov  V.A.,  Romanov  K.V.  Stmctural  analysis  of  Dynamo  zone’s  behaviour  //  Astron.  J.  1993.  Vol.  70. 
P.  880-887. 

7.  Romanov  V.A.,  Romanov  D.V.,  Romanov  K.V.  Magnetic  field’s  emerge  from  Dynamo  zone  up  to  atmosphere  of 
Sun // Astron.  J.  1993.  Vol.  70.  P.  1237-1246. 

8.  Alekseenko  S.V.,  Dudnikova  G.I.,  Romanov  V.A.,  Romanov  D.V.,  Romanov  K.V.  Acoustic  wave  heating  of 
the  Solar  atmosphere  //  Russ.  J.  Eng.  Thermophys.  1998.  Vol.  8.  JV°  1-4.  P.95-108. 

9.  Maoris  C.J.  A  remarkable  eruptive  prominence  on  the  Solar  disk  on  January  29  1968  //  Solar  Phys.  1968.  Vol.  5. 
P.3619-365. 

10.  Altschuler  Martin  D.,  Lilliequist  Carl  G.,  Nakagawa  Yoshinari.  A  Possible  Acceleration  Mechanism  for  a  Solar 
Surge  //  Solar  Phys.  1968.  Vol.  5.  P.  366-376. 


13 


